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1 Abstract 

2 

3 Global ozone analyses, based on assimilation of stratospheric profile and ozone column 

4 measurements, and NOy predictions from the Real-time Air Quality Modeling System 

5 (RAQMS) are used to estimate the ozone and NOy budget over the Continental US 

6 during the July-August 2004 Intercontinental Chemical Transport Experiment- North 

7 America (INTEX-A). Comparison with aircraft, satellite, surface, and ozonesonde 

8 measurements collected during the INTEX-A show that RAQMS captures the main 

9 features of the global and Continental US distribution of tropospheric ozone, carbon 

10 monoxide, and NOy with reasonable fidelity. Assimilation of stratospheric profile and 

1 1 column ozone measurements is shown to have a positive impact on the RAQMS upper 

12 tropospheric/lower stratosphere ozone analyses, particularly during the period when 

13 SAGE III limb scattering measurements were available. Eulerian ozone and NOy budgets 

14 during INTEX-A show that the majority of the Continental US export occurs in the upper 

15 troposphere/lower stratosphere poleward of the tropopause break, a consequence of 

16 convergence of tropospheric and stratospheric air in this region. Continental US 

17 photochemically produced ozone was found to be a minor component of the total ozone 

18 export, which was dominated by stratospheric ozone during INTEX-A. The unusually 

19 low photochemical ozone export is attributed to anomalously cold surface temperatures 

20 during the latter half of the INTEX-A mission, which resulted in net ozone loss during 

21 the first 2 weeks of August. Eulerian NOy budgets are shown to be very consistent with 

22 previously published estimates. The NOy export efficiency was estimated to be 24%, 

23 with NOx+PAN accounting for 54% of the total NOy export during INTEX-A. 

24 

25 
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1 

2 1 . Introduction 

3 One of the key scientific goals of the 2004 phase of the Intercontinental Chemical 

4 Transport Experiment - North America (INTEX-NA) is to quantify and characterize the 

5 inflow and outflow of pollution over North America [Singh, et al. this issue]. The effects 

6 of regional air quality, over the US and elsewhere, on the global atmosphere become 

7 particularly important as world population increases require increases in agricultural 

8 production and continued economic growth leads to increased fossil fuel combustion 

9 [Stevenson et al., 2006]. Combustion leads to anthropogenic emissions of C02, CO, NOx 

10 (NO+N02), S02, and non-methane hydrocarbons (NMHC) as well as particles that can 

1 1 significantly perturb the global atmosphere. In addition to these primary pollutants, 

12 secondary pollutants can have significant impacts on global tropospheric chemistry. In 

13 particular, the abundance and distribution of 03 governs the oxidative capacity of the 

14 troposphere. The global distribution of NOx, which is the critical limiting precursor for 

15 03 production, is highly variable and is dependent on local photochemical loss and 

16 cycling processes involving NOx reservoir species (e.g. PAN, HN03) as well as the 

17 magnitude of various sources which include transport from the stratosphere, natural 

18 emissions (lightning, soils, biomass burning) and anthropogenic emissions (industrial, 

19 aircraft, ships). 

20 This study focuses on estimating ozone and NOy budgets over the Continental US 

2 1 and export to the global atmosphere. The approach we use is an Eulerian budget analysis 

22 as described in Pierce et al. [2003]. This analysis, which focuses on characterization of 

23 the relative contributions of 1) ozone and NOy sources and sinks within the Continental 

24 US domain and 2) regional to global exchange of ozone and NOy, should be interpreted 
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1 in light of the results from a companion manuscript by Al-Saadi et ah, [this issue] which 

2 uses Lagrangian analysis techniques to characterize the influences of the global 

3 atmosphere on the chemical composition of the Continental US. 

4 There is compelling observational and modeling evidence of the link between 

5 Continental US emissions and the global atmosphere. Knapp et al [1998] observed 

6 enhanced ozone (90-130 ppbv) just above the boundary layer over Cape Sable Island, 

7 Nova Scotia during the 1993 North Atlantic Regional Experiment (NARE). Back- 

8 trajectory analysis indicated that these air-masses had origins over the heavily 

9 industrialized N. E. United States. Model studies show episodic but significant remote 

10 influences from North America [Jacob et al., 1993; Wild et al., 1996, Atherton et al, 

1 1 1996, Liang et al., 1998, Li et al. 2004], particularly in the upper troposphere. This 

12 remote influence is driven by export of NOx or PAN, which thermally decomposes to 

13 NOx [Moxim et al., 1996] and leads to further insitu ozone formation [Chameides et al., 

14 1992], Model based estimates of NOy export efficiency suggest that 20-30% of the NOx 

15 emitted from the Continental US is exported to the global atmosphere as NOy 

16 [Kasibhatla et al., 1993; Horowitz et al., 1998, Liang et al., 1998, Li et al., 2004]. 

17 Observational estimates of Continental US NOy export suggest efficiencies ranging from 

18 10-15% [Parrish et al., 2004]. 

19 The preceding discussion illustrates the uncertainties that arise due to complex 

20 interactions between highly heterogeneous surface emissions, local radical chemistry, 

2 1 boundary layer exchange processes, enhancements in background levels of 03 and its 

22 precursors, and long range transport that ultimately determine the links between regional 

23 emissions and the global atmosphere. These links occur across multiple scales in both 
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1 time and space and therefore require a unified approach, utilizing contemporaneous 

2 satellite and insitu observations, as well as model estimates of the chemical state of the 

3 atmosphere. Field missions such as INTEX-A, which use chemical model forecast 

4 guidance to optimize synergy between insitu sampling by airborne platforms and 

5 contemporaneous satellite composition measurements for both satellite validation and 

6 science studies, are an example of this unified approach. However, an “optimized 

7 combination” of satellite, in-situ observations, and model estimates is best accomplished 

8 through chemical data assimilation. Data assimilation provides a physically consistent 

9 representation of the observed atmospheric state and involves blending information from 

10 different sources and different times to yield a best estimate, or “analysis” at a particular 

1 1 time. Models play an important role in data assimilation by providing an estimate, or 

12 “first guess” of the current fields based on previous analyses. The analysis is constructed 

13 by applying an “analysis increment” to the model first guess. The analysis increment is 

14 determined through variational approaches that minimize the differences between the 

15 observation and first guess under constraints that are determined by the relative errors in 

16 the respective fields [Errico, 1999]. 

17 For the current study we utilize ozone analyses (constrained with assimilated 

18 satellite measurements) and NOy predictions from the Real-time Air Quality Modeling 

19 System (RAQMS) [Pierce et ah, 2003] to estimate the ozone and NOy budget over the 

20 Continental US. The manuscript is organized as follows: Section 2 provides an updated 

21 description of the RAQMS, which has undergone significant revisions since Pierce et al., 

22 [2003]. Section 3 focuses on verification of the model 03 analysis as well as CO, NOy, 

23 and 03 P-L predictions based on comparisons with satellite, ozonesonde, airborne, and 
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1 ground based measurements. Section 4 discusses the contributions to ozone and NOy in 

2 the troposphere and lower stratosphere due to stratosphere -troposphere exchange 

3 processes during INTEX-A. The Continental US ozone and NOy budgets during INTEX- 

4 A are presented in Section 5. Section 6 includes a discussion focusing on the 

5 interpretation of the INTEX-A results in light of previous studies. Section 7 provides a 

6 summary and conclusions. 

7 

81. Model Description 

9 The chemical modeling/assimilation tool used in this study is the NASA Langley 

10 Research Center/University of Wisconsin (LaRC/UW) Real-time Air Quality Modeling 

1 1 System (RAQMS). RAQMS is a portable, global- to regional-scale meteorological and 

12 chemical modeling system which has been developed for assimilating remote 

13 observations of atmospheric chemical composition and predicting regional air quality 

14 within any region of the planet Earth [Pierce et al., 2003], This study focuses on the 

15 global modeling/assimilation component of RAQMS. A companion study by Buker et al. 

16 [this issue] utilizes the regional component of RAQMS to investigate stratosphere- 

17 troposphere exchange processes over the Pacific during INTEX-A. The UW hybrid 

18 isentropic coordinate model [Schaack et al., 2004] is the dynamical core for the global 

19 component of RAQMS. Zapotocny et al. [1996, 1997a,b] established that hybrid 

20 isentropic coordinate models simulate processes involving the long-range transport of 

2 1 trace constituents to a higher degree of accuracy than other existing global models. 

22 During INTEX-A RAQMS provided daily 4 day 2x2.5 degree global chemical 

23 forecasts, initialized with ozone analysis based on real-time assimilation of TOMS V8 
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1 ozone column data, to assist in flight planning. The daily assimilation/forecast cycle 

2 consisted of a series of 6 hour online chemical/dynamical forecasts, initialized with 

3 NOAA GFS meteorological analyses at 12Z, 18Z, 00Z, and 06Z. At the end of each 6hr 

4 forecast, the ozone distribution was reinitialized based on the RAQMS TOMS V8 

5 assimilation. After 24 hours of assimilation, a 4-day online chemical/dynamical forecast 

6 was begun. For the current study, we conducted a post mission 1 .4x1 .4 “re-analysis” 

7 from July 01 -August 15, 2004, with meteorological field initialized from the GFS 

8 analyses every 6 hours and including stratospheric ozone profile assimilation in addition 

9 to the TOMS column assimilation. 

10 

1 1 1.1 RAQMS unified stratosphere/troposphere chemistry 

12 The RAQMS unified stratosphere/troposphere chemistry module has been 

13 developed to represent photochemical processes governing ozone formation and 

14 destruction within Earth’s atmosphere from the surface to about 60 km. The chemical 

15 formulation follows a family approach with partitioning based on photochemical 

16 equilibrium approximations. Continuity equations are solved for 55 families and 

17 individual constituents and by determining equilibrium concentrations of 86 separate 

18 species. The standard Ox-HOx-NOx-ClOx-BrOx cycles governing the formation and 

19 destruction of odd oxygen, tropospheric NOx-HOx reactions, oxidation of CH4 and CO 

20 are considered [Pierce et al., 2003], Recent updates include an extended carbon bond 

2 1 scheme for oxidation of non-methane hydrocarbons (NMHC) and explicit treatment of 

22 isoprene oxidation. Photochemical tendencies are calculated with a quasi-steady-state 

23 approximation based on exact solution of the continuity equation. Kinetic rates and 
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1 photolytic quantum yields and absorption cross sections are from Sander et al. [2003] 

2 with the N2+0( l D) quenching rate from Ravishankara et al. [2002]. Photolytic rates are 

3 calculated using the Fastj2 method [Bian et al., 2002]. Stratospheric heterogeneous 

4 reactions on liquid aerosol [Carslaw et al., 1995] and polar stratospheric cloud 

5 [Chipperfield, 1999] surfaces are considered. 

6 The NMHC chemical scheme is based on the lumped-structure approach of the 

7 Carbon Bond -IV mechanism (CB-IV) [Gery et al., 1989] with adjustments necessary for 

8 large-scale (regional or global) application as presented by Zaveri and Peters [1999] 

9 (henceforth called CB-Z). Additional extensions implemented in the LaRC unified 

10 chemistry include an improved isoprene oxidation scheme and the semi-explicit treatment 

11 of propane. The resulting NMHC formulation can be summarized as follows: C2H6 

12 (ethane), C2H4 (ethene) and CH30H oxidation are treated explicitly; C4 and larger 

13 alkanes and C3 and larger alkenes are lumped via a carbon-bond approach as updated in 

14 Zaveri and Peters [1999], with lumped groups for species such as aldehydes, ketones, 

15 peroxides and organic nitrates; C3H8 (propane) is handled semi-explicitly, i.e., C3H8 and 

16 its corresponding peroxy radicals are tracked explicitly (as in Sander et al. [2003] and 

17 Kirchner and Stockwell [1996]) while other oxidation products such as peroxides and 

1 8 aldehydes are lumped into the appropriate species following the carbon-bond approach 

19 for higher alkanes; and isoprene is modeled after the Carter 4-product mechanism as 

20 modified for RADM2 [Carter, 1997]. Aromatic chemistry is not included. Concentrations 

21 of 2 species, acetone and methanol, are currently specified according to climatologies. 

22 Zaveri and Peters [1999] present smog chamber experiments showing that the CB-Z 

23 mechanism is accurate to within 5-20% error under highly polluted conditions. For the 
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1 more remote conditions encountered under regional and global scale modeling, Zaveri 

2 and Peters [1999] compared the CB-Z mechanism to that of RADM2 (Regional Acid 

3 Deposition Model-2 [Stockwell et ah, 1990, 1997]) and found that CB-Z model 

4 calculations were generally within 20% of RADM2 for these multi-day rural simulations, 

5 a significant improvement with respect to the comparison of the original CB-IY with 

6 RADM2 (50-95% differences). The replacement of the CB-Z isoprene oxidation scheme 

7 with a more detailed four-product mechanism allows more accurate representation of 

8 PAN, an important nitrogen reservoir species. The four-product mechanism explicitly 

9 represents the major identified secondary isoprene oxidation products methacrolein 

10 (MACR), methyl vinyl ketone (MVK) and peroxymethacryloyl nitrate (MPAN). Rate 

1 1 constants, products and yields have been updated as reported in the literature (e.g., 

12 Orlando etal. [2002]). 

13 To test the fidelity of the NMHC mechanism within the LaRC unified chemical 

14 module, point calculations using data from the GTE TRACE-P flight campaign were 

15 conducted and compared with results from the LaRC boxmodel [Crawford et ah, 1999]. 

16 The LaRC boxmodel has a NMEIC chemical scheme based on the condensed lumped 

17 molecule mechanism in Lurmann et al. [1986], The two models were run to diurnal 

18 equilibrium with values of 03, CO, NMHCs, temperature, dew point and pressure 

19 constrained to observations. NO was allowed to vary diurnally, but total short-lived 

20 nitrogen (N0+N02+N03+2N205+H0N0+HN04) is held constant such that the 

2 1 predicted NO matched to observations at the time of measurement. Model calculations 

22 are also constrained by H202, CH300H, HN03, PAN and acetone when measurements 

23 are available; otherwise these species are calculated. 
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1 The two models yield nearly identical results when NMHC species are neglected, 

2 thus differences using full chemistries provide a measure of uncertainties introduced by 

3 the choice of NMHC chemical scheme. Median differences between the two models 

4 when NMHC reactions are included are within 10% or less, and the LaRC unified 

5 chemistry is slightly more reactive (lower radical concentrations) compared to the 

6 Lurmann boxmodel mechanism. The one-to-one correspondence between the models is 

7 quite good, with correlation coefficients better than 0.99 for calculated OH, H02 and 

8 CH20. 

9 Climatological emissions of NOx and CO include anthropogenic and natural sources 

10 and are based largely on lxl degree public databases available from GEI A/EDGAR with 

1 1 updates for Asian emissions from Streets et al. [2003] and additional biogenic CO 

12 sources as described by Duncan et al. [2004]. Aircraft NOx emissions are obtained from 

13 the HSRP database [Stolarski et al., 1995]. Lightning NOx emissions are calculated based 

14 on Price et al. [1997] using instantaneous convective cloud heights, and are distributed in 

15 the vertical according to Pickering et al. [1998]. Biomass burning emissions of NOx are 

16 scaled to those of CO. Emissions of NMHC species are generally based on the GEI A 

17 database. Surface sources of N20, CH4, and halocarbons are implicitly assumed by 

18 imposing a constant mixing ratio at the surface appropriate for 1990 [WMO, 1993]. 

19 Surface deposition is computed according to the surface type and drag coefficients, with 

20 the calculation of the deposition rate modeled after Galbally and Roy [1980] and Levy et 

21 al. [1985]. Dry deposition is computed for 03, peroxides, aldehydes, N02, CO and nitric 

22 acid using deposition velocities from Muller and Brasseur [1995]. Online wet removal of 

23 soluble species is based on convective fluxes and precipitation amounts [Liu et al., 2001; 
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1 Park et al., 2004]. Tropospheric heterogeneous loss of N205 is based on zonal averaged 

2 rates from Dentener and Crutzen [1993]. 

3 

4 2.2 RAQMS Ecosystem based wild fire emission inventory 

5 Alaskan and Canadian wild fires had a significant impact on North American 

6 chemical composition during INTEX-A. The RAQMS biomass burning emissions for 

7 Alaska and Canada use an ecosystem based approach developed by Soja et al [2004] to 

8 predict total direct carbon emissions. Soja et al [2004] used a spatially and temporally 

9 explicit model that incorporated a satellite-based (AVHRR) fire database and ecoregion- 

10 specific carbon consumption estimates for three classes of severity to estimate a range of 

1 1 total direct carbon and trace gas emissions from fires in Siberia from 1998 through 2002. 

12 We have extended this algorithm to North America using MODIS thermal anomaly data 

13 to provide area burned estimates [Soja et al., this issue]. Briefly, the approach is as 

14 follows: 

15 (i) Static carbon consumption estimates for low-severity surface fires, medium- 

16 severity fires, and high-severity fires are based on the amount of carbon 

17 contained in individual ecosystems [Olson et al., 1983; Zinke et al., 1986]. 

18 (ii) Daily lxl degree North American area burned estimates are obtained using 

19 MODIS thermal anomaly products. 

20 (iii) Based on carbon consumption and area burned, daily total direct carbon 

2 1 emissions from fire events in Alaska and Canada are estimated for June- 

22 August 2004 assuming all wild-fires were high-severity fires. 
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1 (iv)lxl degree species-specific emission estimates are determined using existing 

2 emission ratios from grassland, temperate and boreal ecosystems from across 

3 North America (Cofer et al., 1996a,b, Vose, 1996). 

4 Assuming high-severity for all Alaskan and Yukon fires must be considered an 

5 upper bound on the actual emissions and results in the release of 70 Tg CO during the 

6 period from June- August 2004, which is a factor of 2 times higher the emissions obtained 

7 by Pfister et al. [2005] during the same period using MOPITT data to constrain inverse 

8 modeling based estimates of the Alaskan/Yukon emissions. In the future, improved 

9 emission estimates that account for changes in fire severity will be obtained using the US 

10 Forest Service Haines Index. The Haines Index is the sum of a stability term and a 

1 1 moisture term. The sum provides an indication of the potential for the rate of spread 

12 (ROS) of a fire on a given day. 

13 

14 2.3 RAQMS Chemical Data Assimilation 

15 Data assimilation provides a statistically robust means of blending model 

16 predictions and observations to provide an optimal estimate of the true state of the 

17 atmosphere. Global assimilation of chemical measurements from polar orbiting satellites 

18 has been shown to improve estimates of the true atmospheric state [Lamarque et al., 

19 1999; Jeuken et al., 1999, Stajner et al., 2004] and is used by RAQMS to provide an 

20 optimal estimate of the global ozone distribution during INTEX-A. RAQMS uses the 

21 statistical digital filter (SDF) analysis system [Stobie 1985, 2000] to perform a univariate 

22 assimilation of stratospheric profile and total column ozone observations. The SDF 

23 formalism is based on optimal interpolation (OI). However, rather than viewing the 
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1 analysis as a minimization problem, SDF treats the analysis as a digital filtering problem 

2 [Oppenheim and Schafer 1975]. In SDF, solving the 01 equations at each grid point is 

3 equivalent to convolving a low-pass digital filter with the observation innovations 

4 (observed value minus first guess value). The spectral response of the filter is determined 

5 by the number of observations used per grid point, the observation spacing, the 

6 observation errors, the first guess and the first guess error correlation model [Stobie 

7 2000], Estimates of the RAQMS forecast error variances are calculated by inflating the 

8 analysis errors (a by-product of the analysis) using the error growth model of Savijarvi 

9 [1995], The quality control employed during the analysis includes a gross check, suspect 

10 identification and a buddy check for suspect observations. 

1 1 Stratospheric (tropopause and above) HALOE, SAGE II, and SAGE III solar 

12 occultation measurements, were assimilated at 6hr (00Z, 06Z, 12Z, 18Z) intervals to 

13 provide constraints on the stratospheric ozone mixing ratios. Assimilation of TOMS V8 

14 cloud cleared total column ozone measurements were used to provide constraints on the 

15 RAQMS total column analysis. The RAQMS column assimilation accounts for the 

16 vertical variation in the retrieval sensitivity by convolving the model first guess ozone 

17 profile with the zonal mean, time averaged sensitivity and the 3D monthly apriori used in 

18 the TOMS V8 retrieval algorithm. Special SAGE III limb scattering measurements 

19 [Rault, 2005; 2006] were taken over North America and the North Atlantic during the 

20 later half of July and mid August. These measurements were also assimilated. Figure 1 

2 1 shows a latitude time series of the frequency of observations used in the RAQMS 

22 assimilation. The symbols indicate the location of the solar occultation and limb 

23 scattering observations while the contours indicate the density of the cloud-cleared total 
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1 column measurements expressed as zonal mean percentages of the total available 

2 observations. During Julian days 183-197 (July 01-14, 2004) SAGE III and HALOE solar 

3 occultation measurements provide profile constraints in the Northern Hemisphere 

4 stratosphere. During Julian days 198-214 (July 15-31, 2004) SAGE III limb scattering 

5 (restricted to North America and North Atlantic sectors) and SAGE II solar occultation 

6 measurements provide additional stratospheric constraints. During Julian days 215-226 

7 (August 01-12, 2004) there are very few stratospheric profile measurements to provide 

8 constraints for the ozone assimilation in the northern hemisphere stratosphere. 

9 

10 3. Model Verification 

11 

12 3.1 Comparison with satellite observations 

13 Figure 2 shows the comparison between RAQMS and cloud cleared MOPITT CO 

14 column for the period from JulyOl-August 15, 2004. The Continental US domain used in 

15 the budget calculations is also shown. To perform this comparison, 6 hourly RAQMS CO 

16 profiles were mapped onto MOPITT observation points and interpolated in time to the 

17 standard MOPITT retrieval levels, then the averaging kernel for each retrieval was used 

18 in conjunction with the MOPITT apriori to determine the “retrieved” RAQMS CO 

19 profile, which was integrated in the vertical using the MOPITT retrieval levels. The 

20 resulting RAQMS “retrieved” and MOPITT retrieved CO columns were binned in lxl 

2 1 degree bins. RAQMS is highly correlated with MOPITT on both regional (continental US 

22 correlation is .92) and global (correlation is .78) scales. RAQMS biases relative to 

23 MOPITT are very small (median biases of less than 0.01 xl0 ls mol/cm 2 , or <1%, globally 
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1 and 0.02xl0 ls mol/cm 2 , or 1.2%, over the Continental US). RAQMS overestimates the 

2 CO column relative to MOPITT over Alaska and Western Canada, where the RAQMS 

3 “retrieved” column is up to a factor of 2 higher than MOPITT. This is consistent with the 

4 factor of 2 higher total wild fire emissions in the RAQMS simulation relative to the 

5 MOPITT constrained emissions used by Pfister et ah, [2005], RAQMS underestimates 

6 the CO column by 50% relative to MOPITT over central Africa. This is because 

7 climatological biomass burning emissions were used in this region. RAQMS also 

8 underestimates the CO column over S.E. Asia and the western Pacific. 

9 Figure 3 shows the comparison between RAQMS and cloud cleared tropospheric 

10 N02 columns retrieved from SCIAMACHY by Martin et al. [in press] for the period 

1 1 from July 01-August 15, 2004. In these comparisons, instantaneous RAQMS N02 

12 profiles were extracted from the model integrations at the SCIAMACHY observation 

13 points to account for the rapid diurnal variation in stratospheric N02. We did not account 

14 for the airmass factor used in the SCIAMACHY tropospheric N02 retrieval. 

15 The AMF calculation of the retrieval uses relative vertical N02 profiles (shape factors) 

16 from GEOS-CHEM. However, little bias is expected in the comparison with RAQMS 

17 since the shape factors are determined largely by the spatial distribution of NOx 

18 emissions [Martin et al., 2002], a distribution which is similar between the two models. 

19 The resulting predicted and measured tropospheric N02 columns were binned in lxl 

20 degree bins. Due to the large dynamic range of the tropospheric column N02 

2 1 measurements the log of the N02 columns are shown. As was found with MOPITT, the 

22 RAQMS spatial distribution is strongly correlated with SCIAMACHY (global correlation 

23 of .65 and continental correlation of .64) but RAQMS tends to underestimate 
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1 tropospheric N02 columns relative to SCIAMACHY both globally (median bias of -.33 

2 xlO 15 mol/cm 2 , or 52%) and regionally (continental US median bias of -.81 xlO 15 

3 mol/cm , or 46%). The systematic low bias over most of the Northern Hemisphere land 

4 masses may be associated with unaccounted for soil emissions [Bertram et al., 2005; 

5 Jaegle et al., 2005]. The relatively coarse resolution of the RAQMS simulation 

6 significantly impacts the inability to accurately describe urban “hot spots”, which also 

7 contributes to the median low biases. 

8 Figure 4 shows the comparison between RAQMS and climatological (1979-2000) 

9 tropospheric ozone determined from TOMS total column and SBUV2 stratospheric 

10 measurements using residual techniques [Fishman and Balok, 1999]. To estimate the 

1 1 climatological July 01 -August 15 mean we have used a 2/3 to 1/3 weighting of the July 

12 and August climatological means. The RAQMS 2004 tropospheric ozone analysis is 

13 generally consistent with climatological expectations both globally (correlation of .73 and 

14 median bias of 1 .09 DU, or 3%) and over the continental US (correlation of .79 and 

15 median bias of 1 .23DU, or 3%) except over Northern Africa and Southern Europe, where 

16 the RAQMS analysis is approximately 10-20 DU higher than climatology. This broad 

17 region of elevated tropospheric ozone column is roughly coincident with the location of 

1 8 the subtropical jet and tropopause break, and was significantly influenced by 

19 stratosphere-troposphere exchange processes during INTEX-A [Al-Saadi et al. this 

20 issue]. 

21 Jing et al. [2004], using contour advection of potential vorticity mapped SAGE II 

22 ozone measurements for 1990, showed enhanced (over 2 Tg/month) isentropic 

23 stratosphere to troposphere ozone transport along the 345K potential temperature surface, 
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1 which is roughly coincident with the mid-latitude tropopause. The enhanced ozone 

2 transport extended in a broad meridional band from the N.E. US, across the Central 

3 Atlantic, and over Northern Africa in a pattern that is remarkably similar to the Atlantic 

4 and European ozone enhancements found in the RAQMS ozone analysis. There is some 

5 evidence of this pattern in the TOMS-SBUV2 climatology, however, the Fishman and 

6 Balok [1999] TOMS-SBUV2 climatology only includes tropospheric ozone residuals that 

7 are less than 75 DU, which would tend to filter out influences of stratosphere to 

8 troposphere exchange processes on the climatological tropospheric ozone. 

9 To determine whether this column ozone enhancement is real or a model artifact 

10 we compared the RAQMS ozone analysis to 8 World Meteorological Organization 

1 1 (WMO) ozonesondes that were launched during the period from July 08 to August 10, 

12 2004 from Santa Cruz, Tenerife, (location shown on Figure 4). Figure 5 shows results 

13 from the statistical analysis of the WMO ozonesonde data. The observed ozone profile 

14 shows persistent ozone enhancements (mean values of 80ppbv) above 500mb. Compared 

15 to the WMO ozonesonde data, the Santa Cruz RAQMS ozone analyses shows mean high 

16 biases of 10% or less below 400mb. However, the high bias in the RAQMS analysis 

17 increases to nearly 35% at the tropopause (near 150mb), suggesting that while there is 

18 clearly upper tropospheric ozone enhancement at Santa Cniz it is overestimated in the 

19 RAQMS ozone analysis. A 35% overestimate in upper tropospheric ozone has a 

20 relatively small effect on the tropospheric ozone column. The mean observed and 

21 analyzed ozone columns below 150mb were 43.9DU and 48.4DU, respectively. The 

22 resulting 4.5DU overestimate is approximately 10% of the observed column. 

23 
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1 3.2 Comparison with IONS ozonesonde and EPA AIRNOW networks 

2 The INTEX Ozonesonde Network Study (IONS) [Thompson et al, this issue] 

3 provided multiple daily ozonesonde launches during INTEX-A. The unprecedented 

4 duration (July 01-Aug 14), frequency (daily), and density (up to 12 ozonesonde stations) 

5 makes this data set extremely useful for verification of the RAQMS ozone assimilation as 

6 well as science studies [Thompson et al, 2006]. Figure 6 shows composite timeseries of 

7 the IONS ozonesonde data along with comparisons with the RAQMS ozone analysis. The 

8 IONS composite was obtained by binning all the daily ozonesondes in log-pressure bins 

9 (10 bins/decade in pressure). The RAQMS composite was obtained in the same manner 

10 after mapping the RAQMS ozone analysis to the full resolution ozonesonde. Mean errors 

1 1 were estimated by averaging the point-by-point errors for all profiles within each pressure 

12 bin. RMS errors were estimated in the same way after removing the mean bias within 

13 each pressure level. The daily mean pressure of the thermal tropopause on the northern 

14 boundary of the Continental US budget domain, and 380K potential temperature surface 

15 on the southern boundary of the budget domain are also indicated. These surfaces define 

16 the maximum vertical extent of the middle world [Holton, 1999] over the Continental 

17 US. The middle world is a region in the lower stratosphere that is strongly coupled to the 

1 8 subtropical upper troposphere due to quasi-horizontal, isentropic exchange near the 

19 subtropical jet. This quasi-horizontal exchange occurs because of the sharp meridonal 

20 gradient in the tropopause altitude near the subtropical jet. The IONS composite shows 

2 1 significant day to day variability in this region, as evidenced by the altitude of the 

22 250ppbv ozone mixing ratio (which is roughly coincident with the thermal tropopause) 

23 prior to mid July (Julian day 200) then there is an extended period with less variability, 
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1 followed by renewed variability during early August (Julian day 215). The RAQMS 

2 ozone analysis does a good job in reproducing this composite behavior, which is largely 

3 driven by upper tropospheric planetary wave activity. 

4 Upper tropospheric ozone mixing ratios of 100 ppbv or more are most likely 

5 stratospherically influenced airmasses. The IONS composite shows significant variability 

6 in the frequency of stratospherically influenced tropospheric ozone measurements, 

7 defined here as composite ozone mixing ratios greater than 100 ppbv but observed below 

8 northern tropopause of the Continental US domain. In the IONS composite, signatures of 

9 stratospherically influenced ozone extend down to 500mb during the middle 3 weeks of 

10 INTEX-A. There are relatively fewer observations of stratospheric influenced air in the 

1 1 troposphere during the first 2 weeks of July (Julian days 183-198) and second week of 

12 August (Julian days 227-234). This is consistent with Avery et. ah, [this issue] who find 

13 evidence for interleaving and mixing of stratospherically influenced and polluted 

14 tropospheric air in the vicinity of the subtropical jet for INTEX-A flights in late July and 

15 early August. 

16 The observed variations in the depth of the stratospherically influenced air in the 

17 upper troposphere are not as pronounced in the RAQMS composite, leading to 

18 predominately positive analysis errors in the upper troposphere during the beginning and 

19 end of INTEX-A and predominately negative analysis errors in the upper troposphere 

20 during the middle of INTEX-A. The positive analysis errors at the beginning an end of 

2 1 the mission frequently reach 40-50% and extend down to as far as 600mb while the 

22 negative analysis errors are typically only 20-30%. During the last two weeks of July and 

23 first week of August the frequency and vertical extent of large positive analysis errors are 
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1 significantly reduced , with the majority of the analysis errors between +/- 20% during 

2 the middle part of INTEX-A. 

3 The first part of this period of relatively low analysis errors (July 15-3 1) 

4 corresponds to the period when daily SAGE III limb scattering measurements, made over 

5 a wide latitude band over the Continental US, were assimilated. The reduction in the 

6 extent and frequency of significant high biases in the analysis during this period indicates 

7 that the assimilation of SAGE III limb scattering measurements had a positive impact on 

8 the RAQMS ozone analysis, particularly in the upper troposphere/lower stratosphere. The 

9 fact that these improvements persist for at least 5 days after the limb scattering 

10 assimilation stops indicates that the system has memory of the measurements, which has 

1 1 significant implications for air quality forecasting. The RMS analysis errors over the 

12 Continental US are largest within the middle world and are typically on the order of 40- 

13 60%, although occasionally RMS errors reach up to 80-100%. These errors are most 

14 likely associated with errors in the vertical placement of stratospheric and subtropical 

15 ozone lamina that were frequently observed during INTEX-B [Thompson et al., 2006], 

16 There are occasional RMS errors of up to 40% near the surface, but these appear to occur 

17 during periods of relatively low boundary layer ozone events and are not likely to be 

18 significant. 

19 Figure 7 summarizes the RAQMS/IONS comparison with a time averaged 

20 comparison between the RAQMS ozone analysis and all ozone profiles during INTEX-A. 

21 In addition to mean and rms errors we also assess the ability of the RAQMS analysis to 

22 capture the observed variability, defined here as the overall temporal and site to site 

23 variability at a given pressure level. Above lOOmb, the mean analysis biases are on the 
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1 order of 10% with rms errors of less than 20%. The upper troposphere/lower stratosphere 

2 shows mean high biases of near 20%, extending from 100mb-300mb, below 300mb, the 

3 RAQMS shows low biases of less than 10%. The RAQMS analysis captures the majority 

4 of the observed variance enhancement in the lower stratosphere and continental boundary 

5 layer but also shows large RMS errors (near 50% at 200mb). These RMS errors are due 

6 to vertical displacement of fdaments of high and low ozone associated with stratosphere - 

7 troposphere exchange processes. 

8 Figure 8 shows comparisons between RAQMS surface ozone and ozone 

9 measurements from the EPA AIRNOW network [Wayland, 2002]. These maps show 

10 mean statistics based on timeseries analysis for the individual AIRNOW stations. The 

1 1 median temporal correlation between the 6 hourly RAQMS prediction and coincident 

12 lhour AIRNOW measurements is .72 (not shown), and largely reflects the diurnal cycle 

13 in surface ozone. To assess the ability of the RAQMS ozone analysis to capture daily 

14 variations in surface ozone we consider correlations between the diurnally averaged 

15 RAQMS analysis and coincident AIRNOW measurements. The median correlation 

16 between diurnally averaged RAQMS and AIRNOW data is .59 with lowest correlations 

17 in the diurnally averaged ozone found over West Virginia, Southern California, and the 

18 western mountain states. These low correlations are most likely associated with 

19 unresolved local variations in dry deposition due to local variations in topography and 

20 emissions. 

21 To assess the ability of the RAQMS ozone analysis to capture daytime 

22 photochemistry we compiled station-by- station mean biases at 18Z, which is mid-day 

23 over much of the central and eastern US. 18Z mean biases are generally positive with a 
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1 median value of 15.5 ppbv. 18Z mean biases are largest within the Mississippi and Ohio 

2 River valleys. The daytime mean biases may associated with overestimates in surface 

3 ozone production, overestimates in boundary layer 03 entrainment, assimilation of 

4 TOMS column ozone (which occurs during the 18Z assimilation cycle over North 

5 America), or overestimates in the initial (morning) surface ozone. Overestimates in 

6 surface ozone production would suggest excess NOx, however column N02 is actually 

7 underestimated (particularly in urban areas) based on the RAQMS/SCIAMACHY 

8 comparison. Entrainment is also unlikely to account for the 1 8Z bias since boundary layer 

9 03 is actually underestimated (in the mean) based on the RAQMS/IONS comparison. 

10 This leaves overestimates in the initial (morning) surface ozone or TOMS column ozone 

1 1 assimilation as the most likely reasons for the mean daytime bias. Nighttime (OOZ) 

12 biases show a similar pattern with median values of 18.1 ppbv. The OOZ biases are 

13 frequently associated with underestimates in nighttime titration of ozone, which is a near- 

14 surface phenomena that is not accurately captured within RAQMS. 

15 

16 3.3 Comparison with DC8 insitu measurements 

17 Figure 9 shows comparisons with insitu CO, 03, N02, PANS, and HN03 data 

18 obtained by instruments onboard the NASA DC8 during all flights during INTEX-A as 

19 well as observationally constrained photochemical steady state calculations from the 

20 LaRC boxmodel [Crawford et ah, this issue]. These comparisons were made by 

21 interpolating the RAQMS chemical fields onto the DC8 flight track and sampling the 

22 model at the frequency of the insitu measurements, then binning the modeled and 

23 measured values into 50mb pressure bins. The median (vertical profile), 50 th (bar) and 
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1 90 th (whisker) percentiles of the modeled and observed distributions within each pressure 

2 bin are shown. The modeled CO is 5-10 ppbv lower than observed except at 950mb 

3 where the model is approximately 20 ppbv higher than insitu measurements. Predicted 

4 and observed column amounts, obtained by integrating the median number densities 

5 between 1000 and 250 mb, are 1.8xl0 18 and 1.82xl0 18 mol/cm 2 , respectively. The 1-2% 

6 agreement between the predicted and insitu median column amounts is consistent with 

7 the comparison between the RAQMS and MOPITT CO columns, which showed median 

8 biases of 1.2% over the continental US. The modeled 03 is within lOppbv below 400mb 

9 and approximately 20ppbv higher than the observations above 400mb. The column 

10 densities obtained from integration of the median analyzed and insitu ozone profdes are 

1 1 46.6 and 42.2 DU, respectively, resulting in a 10% error in the median ozone column 

12 estimated from the aircraft measurements. The 20-25% differences between the RAQMS 

13 ozone analysis and insitu mixing ratios in the upper troposphere are consistent with the 

14 comparison between RAQMS and the IONS ozonesonde data, suggesting that the 

15 analyzed tropospheric 03 column is within 10% of the actual column over the continental 

16 US. The modeled N02 is approximately 50ppbv low (factor of 2) relative to the 

17 measurements at 300mb and larger than observed below 800mb. The predicted median 

18 N02 mixing ratio is a factor of 2 larger than observed at 950mb. Column N02 densities, 

19 based on integration of the predicted and observed median profiles, are 1.74xl0 15 and 

20 1.34x10 mol/cm", respectively, resulting in a 30% bias in the median N02 column 

2 1 estimated from the aircraft measurements. The high bias in median N02 column density, 

22 relative to the insitu data, is not consistent with the 46% negative median bias found 

23 relative to SCIAMACHY tropospheric N02 column densities. The differences between 
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1 the results of the SCIAMACHY and insitu verification studies could only arise from low 

2 biases in the RAQMS N02 mixing ratios below the 1000ft minimum altitude of the DC8. 

3 Such biases are very likely within urban boundary layers (which were not sampled by the 

4 DC8) and could also result from underestimates in large soil NOx emissions associated 

5 with fertilizer application and subsequent precipitation [Bertram, et ah, 2005, Juegle, et 

6 ah, 2005] in the western US. Total PANs (peroxynitrates) are within 50pptv of the 

7 observed mixing ratios except for near 800mb where the modeled Total PAN 

8 (PAN+HN04) is low by 75pptv. HN03 is low by 200 - 400 pptv below 600mb and high 

9 by 150 pptv at 300mb. The modeled ozone P-L is in very good agreement with 

10 observationally constrained photochemical steady state calculations except at 300mb 

1 1 where the modeled P-L (1 ppbv/day) is low by a factor of 2. This underestimate in ozone 

12 P-L is consistent with factor of 2 underestimate of N02 at this altitude. The fact that the 

13 model ozone overestimates are associated with overestimates in HN03 and 

14 underestimates in P-L, N02 and TOTPANS, suggests that the model overestimates 

15 stratospheric influences (high 03, HN03), underestimates convective influences (P-L, 

16 N02, PAN), or both in the upper troposphere. 

17 

18 4. Global and Continental US estimates of Ozone and NOy STE 

19 The preceding discussion highlights the important role that stratospheric- 

20 tropospheric exchange (STE) plays in determining the distribution of ozone and NOy 

2 1 (primarily HN03) in the upper troposphere during INTEX-A. In this section, we follow 

22 the discrete approach outlined by Pierce et al, [2003] to estimate the contribution of STE 

23 to the global distribution of upper tropospheric ozone during INTEX-A. In Pierce et al., 
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1 [2003] the regional component of RAQMS was used to determine discrete, cross- 

2 tropopause ozone fluxes over S.E. Asia during the NASA TRACE-P mission. The 

3 discrete cross tropopause flux was estimated by computing instantaneous horizontal and 

4 vertical fluxes out of tropospheric grid-boxes that were adjacent to the model tropopause, 

5 using the WMO thermal tropopause definition. Here we apply the same approach using 

6 the global component (UW-Hybrid dynamical core) of RAQMS. The UW-Hybrid model 

7 is formulated in hybrid isentropic-eta coordinates and consequently grid boxes are 

8 defined in the vertical by potential temperature surfaces above 345K, or roughly the mid- 

9 latitude tropopause. This hybrid isentropic-eta formulation of the UW-Hybrid dynamical 

10 core allows us to explicitly compute the isentropic (quasi-horizontal) exchange of 

1 1 stratospheric and tropospheric air across the tropopause break, which extends from 

12 roughly 345K to 380K in potential temperature. This region of the lower stratosphere, 

13 bounded by the tropopause and the 380K potential temperature surface, is coupled to the 

14 upper tropical troposphere through isentropic exchange of mass, momentum, and trace 

15 gases. 

16 Figure 10 shows the zonally averaged cross tropopause ozone and NOy fluxes 

17 during the period from July 01- August 15, 2004. The ozone and NOy fluxes across the 

18 380K potential temperature are also shown. NOy fluxes are determined by computing the 

19 instantaneous 6 hourly fluxes of individual components of NOy 

20 (NO+N02+N03+HN03+HN04+2*N205+ClN03+PAN+organic nitrates) and then 

2 1 adding them together. The RAQMS first guess odd oxygen (Ox) is used to compute the 

22 6hourly ozone fluxes to assure dynamical consistency between the Ox and forecasted 

23 winds. The time averaged cross tropopause ozone and NOy fluxes is determined by 
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1 averaging 6 hourly calculations of horizontal (isentropic) and vertical (diabatic) fluxes, 

2 and movement (in altitude) of the tropopause. Since the 380K surface is a model level, 

3 the horizontal velocities at 380K are parallel to the 380K surface and the time averaged 

4 ozone and NOy flux at 380K only includes vertical fluxes plus movement (in altitude) of 

5 the 3 80K surface. 

6 In the tropics, the ozone and NOy fluxes are upward (positive) with net transport 

7 through the tropopause and 380K potential temperature surface into the lower 

8 stratosphere. This transport is driven by radiative heating and upward diabatic vertical 

9 motion. In polar regions, the diabatic transport of ozone and NOy across the 380K 

10 potential temperature surface is downward and driven by radiative cooling. The net 

1 1 transport of ozone and NOy across the polar tropopause is also downward, but it occurs 

12 through the combined effects of downward diabatic motion and secular changes 

13 (increases) in the altitude of the polar tropopause, which compensate for net horizontal 

14 (isentropic) fluxes of ozone and NOy into the middle world. At the tropopause break, 

15 horizontal (isentropic) ozone and NOy fluxes from the troposphere into the stratosphere 

16 dominate. 

17 These results are consistent with mass flux estimates by Schoeberl [2004] who 

18 used explicit calculations of the diabatic fluxes through the 380K potential temperature 

19 surface and tropopause, along with mass tendencies within the middle world, to estimate 

20 adiabatic (isentropic) exchange through the tropopause. He found net adiabatic fluxes 

2 1 across the tropopause were positive (into the middle world) throughout the year. 

22 However, because the adiabatic term was obtained as a residual, he was not able to 

23 determine the latitudinal distribution of the adiabatic cross tropopause flux. The RAQMS 
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1 analysis shows that the adiabatic flux is largest at the tropopause break. This net flux 

2 (from above and below) of trace gases into the middle world on the poleward side of the 

3 tropopause break introduces the possibility for accumulation of ozone and NOy with both 

4 stratospheric and tropospheric origins within this region. 

5 Figure 1 1 shows the zonally averaged cross tropopause ozone and NOy fluxes 

6 during the period from July 01- August 15, 2004 for the Continental US budget domain. 

7 The regional ozone and NOy fluxes across the 380K potential temperature surface are 

8 also shown. The cross tropopause fluxes of ozone and NOy over the Continental US are 

9 dominated by quasi-horizontal (isentropic) transport into the middle world. The both 

10 ozone and NOy show peak troposphere to stratosphere ozone fluxes at 45oN with the 

1 1 cross tropopause NOy flux being dominated by transport of HN03. The fluxes of ozone 

12 and NOy across the 380K potential temperature surface are largest on the northern 

13 boundary of the budget domain and are dominated by downward (diabatic) transport into 

14 the middle world. These results are consistent with the global flux estimates for this 

15 latitude band and suggest that there should be an accumulation of ozone and NOy, some 

16 of tropospheric and some of stratospheric origin, in the lower stratosphere over the 

17 Continental US during INTEX-A. The net upward flux of ozone and NOy at the 

18 tropopause within the Continental US budget domain suggests that the stratospherically 

19 influenced tropospheric air observed in the IONS composite must have entered the 

20 troposphere poleward of the Continental US domain, where net cross tropopause fluxes 

2 1 are downward. 

22 

23 5. Ozone and NOy budgets over the Continental US 


27 



1 5.1 Ozone Budget analysis 

2 Ozone assimilation reduces errors in the budget analysis by providing an 

3 improved estimate of ozone within the Continental US budget volume. However, 

4 assimilation introduces non-physical changes in ozone that must be isolated from the 

5 physical and chemical processes accounted for in the RAQMS simulation. To isolate the 

6 influences of assimilation in the budget calculations we use the RAQMS first guess 

7 ozone distributions to compute the lateral and diabatic fluxes. Since the first guess ozone 

8 has been advected for the previous 6 hours with forecasted wind fields and experienced 

9 the effects of the predicted photochemistry, the first guess is dynamically and chemically 

10 consistent. The ozone analysis increment is treated as a separate (although non-physical) 

1 1 budget term. This approach allows us to isolate the effects of assimilation in the budget 

12 calculations. 

13 Figure 12 shows the time averaged zonal mean distribution of Continental US 

14 ozone, net P-L, convective mixing tendencies, and the mean absolute value of the 

15 assimilation increment. The mean location of the middle world, bounded by the thermal 

16 tropopause and 380K potential temperature surface, is also indicated. Mean ozone mixing 

17 ratios are above 80ppbv in the northern upper troposphere and range from 200-350 in the 

18 middle world. Lower tropospheric ozone mixing ratios are less than 55 ppbv, with no 

19 clear indication of a surface enhancement. Time averaged upper tropospheric net ozone 

20 production (P-L) reaches 4.0 ppbv/day at 10 km and 30N, and shows net photochemical 

2 1 destruction below 7 km in the southern portion of the domain. Net photochemical 

22 production reaches lOppbv/day in the Continental US boundary layer. Deep convection 

23 in the southern part of the domain leads to upper level detrainment of low ozone mixing 
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1 ratios within convective updrafts, resulting in upper tropospheric ozone reductions of 

2 nearly 7ppbv/day. This convective ozone sink is localized near the region of largest net 

3 photochemical production, suggesting lighting NOx emissions and convective transport 

4 of boundary layer NOx emissions play an important role in the upper tropospheric ozone 

5 production. Lateral detrainment of higher ozone mixing ratios during deep convection 

6 leads to mid-tropospheric increases in ozone at a rate of 5 ppbv/day. Shallow convection 

7 near 40N results in entrainment of higher ozone mixing ratios associated with ozone 

8 production within the continental boundary layer. This low level entrainment and 

9 subsequent convective lofting leads to localized convectively induced ozone sinks of up 

10 to 2.0 ppbv/day in this region. The RMS effects of the ozone assimilation are less than 

11 1% over much of the troposphere, with localized regions below 2km in the northern part 

12 of the budget domain showing upwards of 2% RMS changes. Relatively uniform 

13 assimilation increments of 2% or more are found just above the tropopause and are a 

14 result of the assimilation of the solar occultation and limb scattering measurements. 

15 Figure 13 shows the time averaged ozone number densities for each of the lateral 

16 boundaries of the Continental US domain. The time averaged tropopause and 380K 

17 potential altitudes are also shown. On the western boundary the ozone number densities 

1 8 are very low below 2km, reflecting the influence of clean maritime air. Mid-tropospheric 

19 ozone enhancements are found along each domain boundaries. Within the middle world, 

20 there are local maxima in ozone number densities near the northern edges of both the 

2 1 western and eastern boundaries, as well as over the western half of the northern 

22 boundary. These local maxima are evidence for the accumulation of ozone within the 

23 middle world and are a result of the global scale vertical and horizontal flux convergence 
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1 poleward of the tropopause break discussed in Section 4. The largest local enhancements 

2 in ozone number densities are found in the middle world along the western portion of the 

3 northern boundary. These ozone enhancements are associated with an upper tropospheric 

4 trough pattern which persisted during much of INTEX- A. The signature of this upper 

5 tropospheric trough is evident in the lower time averaged tropopause altitudes along the 

6 western portion of the northern boundary. 

7 Figure 14 shows the time averaged ozone fluxes, in mol/cm 2 /sec, for each of the 

8 lateral boundaries of the Continental US domain. Negative values denote fluxes into the 

9 Continental US while positive values denote export out of the Continental US. The 

10 largest ozone fluxes (both into and out of the Continental US) occur along the eastern and 

1 1 western boundaries and arise due to the prevailing westerly winds along the northern 

12 portion of the Continental US. These ozone fluxes maximize in the middle world, and are 

13 coincident with local maxima in ozone number densities shown in Figure 13. There is a 

14 reversal of the fluxes on the eastern and western boundaries over the southern US 

15 associated with the prevailing stratospheric easterlies near 20km. The ozone export along 

16 the eastern boundary is significantly larger than ozone import along the western 

17 boundary, leading to net ozone export out of the Continental US. The alternating pattern 

18 of middle world ozone fluxes along on the northern boundary is a consequence of 

19 meridional transport within the upper level trough, with net import on the western flank 

20 and export on the eastern flank. The ozone fluxes along the southern boundary are very 

2 1 small. 

22 The large time averaged lateral fluxes, coupled with flux convergence associated 

23 with upward STE and downward diabatic transport into the middle world, suggest that 
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1 neglecting this region in the US ozone budget could lead to significant underestimates in 

2 the actual export during INTEX-A. Consequently, in the subsequent budget analysis, we 

3 consider the 380K potential temperature surface to be the top of the budget domain. 

4 Figure 15 shows the time series of the accumulated changes in Continental US ozone 

5 associated with ozone production, 380K diabatic fluxes, lateral fluxes, ozone 

6 assimilation, and dry deposition. The actual and computed accumulation, determined 

7 from the sum of the individual budget terms, is also shown. The initial Continental US 

8 ozone burden below 380K was 15.6 Tg. Variations in the total ozone below 380K are on 

9 the order of 2 Tg with small (<1 Tg) net changes in Continental US ozone during 

10 INTEX-A. The close agreement between the actual and computed accumulation during 

1 1 INTEX-A indicates low accumulative errors in the budget calculation. 

12 Net photochemical production is the dominant source of changes in Continental 

13 US ozone during INTEX-A, with accumulated insitu ozone production of 7.63 Tg. 

14 However, photochemical production slows down significantly after Julian Day 198 (July 

15 15), as reflected in the much slower accumulation due to ozone photochemistry during 

16 the later half of July, and becomes negative in August, as reflected in the decline in ozone 

17 accumulation due to Net P-L after Julian day 220. The reductions in photochemical ozone 

18 production are attributed to the anomalously cold weather pattern during these periods. 

19 August 2004 was the 7 th coldest on record with cold outbreaks occurring during July 26- 

20 30 and August 10-16 [Fuelberg et al, this issue]. These cold air outbreaks are associated 

2 1 with increased surface winds which lead to efficient boundary layer ventilation and 

22 reduced accumulation of ozone precursors. US EPA air quality statistics for 92 major 

23 metropolitan areas in the continental US show that 2004 had the fewest days with ozone 
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1 Air Quality Indexes over 100 (corresponding to 8 hour average ozone mixing ratios 

2 greater than 85 ppbv) during the last 15 years (1990-2004). The reduction was highly 

3 significant. When all US metropolitan areas are considered the number of ozone 

4 AQI>100 days in 2004 was 66% less than the 15 year median value. If we exclude 

5 California, Dallas, and Houston the number of ozone AQI>100 days in 2004 was 82% 

6 less than the 15 year median. 

7 The rate of ozone loss due to dry deposition remains nearly constant throughout 

8 the INTEX-A period and is the dominate sink of ozone within the budget volume. 

9 Accumulated losses due to dry deposition (7.39 Tg) nearly balance net photochemical 

10 production over the Continental US during INTEX-A. Lateral ozone fluxes result in 

1 1 accumulative reductions (net export) of 3 Tg during INTEX-A, however, most of this 

12 export occurs prior to Julian Day 201 (July 18 th ). After July 18 th there were a series of 

13 anomalously deep upper level troughs over the US [Fuelberg, et al, this issue] which 

14 contributed to the significant week to week variability in the lateral fluxes during the 

15 latter half of July and first half of August. Diabatic fluxes across the 380K surface result 

16 in the import of 3 Tg of stratospheric ozone to the Continental US domain during the 

17 INTEX-A time period. 

1 8 There is a significant reduction in the rate of accumulation of stratospheric ozone 

19 within the Continental US domain after Julian day 197 (July 15 th ) corresponding to the 

20 beginning of the assimilation of SAGE III limb scattering measurements. As was shown 

21 in Section 3.2, assimilation of SAGE limb scattering measurements had a positive impact 

22 on the RAQMS vs IONS ozonesonde statistics by reducing high biases in the region of 

23 the upper troposphere with significant stratospheric influences. Assimilation increments 
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1 during the SAGE limb scattering period (July 15 th - July 3 1 st ) result in a net loss of 3 Tg of 

2 ozone over the continental US, which is comparable in magnitude to the total ozone 

3 export during INTEX-A. After July 3 1 st the assimilation of the SAGE limb scattering 

4 measurements stops and assimilation of TOMS column ozone only results in systematic 

5 increases in ozone within the budget volume until the SAGE limb scattering assimilation 

6 is resumed on August 12 th . The assimilation of SAGE III limb scattering measurements 

7 also impacts the estimates of lateral fluxes, which maximize in the middle world. This is 

8 reflected in the anti-correlation between the accumulated effects of assimilation and 

9 lateral fluxes after July 15 th . Whether changes in the lateral fluxes are associated with the 

10 inclusion of SAGE III limb scattering data in the assimilation or changes in the upper 

1 1 tropospheric circulation during the latter part of INTEX-A is a difficult question to 

12 answer. Future budget studies could address this question by conducting budget analysis 

13 with and without assimilation. 

14 The majority of the export from the Continental US domain occurs in the middle 

15 world, consequently, the net import of approximately 3 Tg of stratospheric ozone across 

16 the 380K potential temperature surface is likely to account for the majority of the 3 Tg of 

17 ozone that is exported during INTEX-A. To obtain an estimate of the export of ozone that 

18 was photochemically produced within the Continental US domain during INTEX-A we 

19 need to remove the stratospheric contribution from the accumulated lateral export. This 

20 results in 9.4e9 g of ozone photochemically produced over the US and exported during 

2 1 INTEX-A, which is a negligible fraction of the total export. 

22 

23 5.2 NOy Budget analysis 
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1 The export of total reactive nitrogen (NOy) from the continental US is equally as 

2 important as the export of ozone, since availability of nitrogen oxides (NO+NO?) 

3 determine subsequent ozone production [Chameides et ah, 1992]. In this section we 

4 discuss results from Eulerian budget calculations focusing on NOy. Figure 16 shows the 

5 time averaged zonal mean distribution of Continental US NOy, lighting NOx production, 

6 convective exchange of NOy, and NOy wet deposition. The zonal mean surface NOy is 

7 over 4 ppbv and is dominated by localized NOx enhancements due to emissions and 

8 HN03. There is a pronounced tongue of elevated NOy extending down from the mid- 

9 latitude tropopause that has significant stratospheric influences. Since the cross 

10 tropopause NOy flux is from the troposphere to the stratosphere within the Continental 

1 1 US budget domain, these NOy enhancements must arise due to STE outside of the 

12 Continental US. Al-Saadi et al. [this issue] show much of the stratospherically influenced 

13 air within the continental US had it’s origins over the central Pacific and S. E. Asia. 

14 Cloud top detrainment of lightning NOx emissions in the southern portion of the domain 

15 results in NOy production of 0.25 ppbv/day in the upper troposphere (7- 10km), with 

16 nearly equal amounts below 2km associated with outflow from convective downdrafts 

17 [Pickering, 1998]. Convective mixing entrains continental boundary layer NOy at a rate 

18 of 2 ppbv/day where it is either immediately rained out (for highly soluble species such 

19 as HN03) or convectively lofted (for less soluble species such as PAN) and deposited 

20 between 5 and 10 km. In contrast to ozone, convective exchange increases free 

2 1 tropospheric NOy mixing ratios at a rate of 0. 1 5 ppbv/day. This is a consequence of the 

22 different vertical gradients in NOy and ozone below 10 km. 
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1 Figure 17 shows the time averaged NOy number densities for each of the lateral 

2 boundaries of the Continental US domain. The distribution of middle world NOy and 

3 ozone number densities (Figure 13) are very similar with local maxima on the northern 

4 edges of the western and eastern boundaries and on the eastern edge of the northern 

5 boundary. These local maxima are primarily HN03, and result from net flux converge 

6 within this region. However, in the troposphere there are significant differences between 

7 the ozone and NOy number densities on the lateral boundaries. The largest NOy number 

8 densities are found below 5km on the northern boundary and below 2 km on the eastern 

9 boundary. The enhancements in NOy on the northern boundary are primarily HN03 and 

10 PAN and are due to transport from the Alaskan wild fires [Al-Saadi et al, this issue]. The 

1 1 large local enhancements in NOy number densities below 2km along the eastern 

12 boundary are primarily due to HN03 as are the low level enhancements in NOy on the 

13 western boundary. Mid tropospheric enhancements in NOy number densities along the 

14 eastern and southern boundaries are primarily due to PAN. Figure 18 shows the time 

15 averaged NOy fluxes, in mol/cm7sec, for each of the lateral boundaries of the 

16 Continental US domain. As with the number densities, the distribution of NOy fluxes in 

17 the middle world are similar to the ozone fluxes (Figure 14) and are dominated by fluxes 

18 of HN03. In the troposphere, there is significant NOy import (negative fluxes) on the 

19 northern boundary near 5km. These fluxes are primarily due to transport of HN03 and 

20 PAN from the Alaskan wildfires. The NOy export (positive fluxes) on the northern part 

21 of the eastern boundary extends well into the troposphere. This is due to export of PAN, 

22 which maximizes near 7km along the northern portion of the eastern boundary. The 

23 localized export of NOy below 2 km at 45N is primarily composed of HN03. 


35 



1 While there is significant complexity in the way that NOy species are partitioned 

2 among the various regions of import and export the evolution of the accumulated 

3 changes in Continental NOy is actually quite simple due to the fact that it’s primary 

4 source is surface emissions, which are held constant throughout the simulation. Figure 19 

5 shows the time series of the accumulated changes in Continental US NOy (expressed in 

6 Tg of nitrogen) due to sources (industrial plus aircraft and soil emissions, lightning 

7 emissions), sinks (wet and dry deposition) and transport (380K diabiatic fluxes and 

8 lateral fluxes). The actual and computed accumulation, determined from the sum of the 

9 individual budget terms, is also shown. The initial Continental US NOy burden below 

10 380K was 0.05 Tg. The actual and computed NOy accumulation over the Continential 

1 1 US were very small during 1NTEX-A, as are the 380K diabatic fluxes. The Continental 

12 US NOy budget shows accumulated NOy emissions of 0.94 Tg nitrogen (with less than 

13 20% due to lightning NOx production) and accumulated depositional loss of 0.69 Tg 

14 nitrogen (.47 Tg wet, .22 dry), resulting in a net export of 0.23 Tg of nitrogen and an 

15 export efficiency of 24%. 

16 

17 6. Discussion 

18 Liang et al. [1998] (here after referred to as L98) used sensitivity experiments 

19 based on differences between two continental-scale photochemical model simulations 

20 (one with and one without US NOx emissions) to estimate seasonally averaged fluxes of 

2 1 ozone and NOy. The summer season (JJA) ozone export from the Continental US 

22 boundary layer in the standard (with US NOx) simulation was 1.8 Gmol/day, while 

23 difference between the standard simulation and one without US NOx emissions, referred 
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1 to as “pollution ozone” was 6.5 Gmol/day. Li et al. [2004] (here after referred to as L04) 

2 used the GEOS-CHEM model sensitivity experiments to estimate US ozone export out of 

3 the Continental US boundary layer during September 1997 and found “pollution ozone” 

4 export of 5 Gmol/day, consistent with the Fall (SON) estimates by L1998. Furthermore, 

5 L04 showed that nearly 70% of the ozone production associated with Continental US 

6 NOy export out of the boundary layer occurs directly over North America, referred to as 

7 “near field ozone production”, and would therefore be included in the RAQMS INTEX-A 

8 Continental US budget calculations presented here. 

9 The L98 horizontal Continental US domain size was similar to the current 

10 INTEX-A budget domain. Applying the JJA L98 standard simulation export rate over the 

1 1 46 day INTEX-A budget period would result in 3.97 Tg of ozone exported through the 

12 continental US boundary layer, comparable to our estimates of net export, however, as 

13 discussed earlier, the RAQMS INTEX-A ozone export is of stratospheric origin. 

14 Applying the JJA L98 “pollution ozone” export rate over the 46 day INTEX-A budget 

15 period would result in 14.35 Tg of US ozone exported through the continental boundary 

16 layer, which is significantly larger than our estimates of net export of photochemically 

17 produced ozone from the Continental US budget domain during INTEX-A. 

18 Direct comparisons of the INTEX-A Continental US photochemical ozone export 

19 and the L98 and L04 “pollution ozone” estimates are not appropriate since the “pollution 

20 ozone” reflects the fact that without emissions, the Continental US would be a strong sink 

21 of ozone due to photochemical losses and dry deposition near the surface. However, 

22 because of the large discrepancies between the L98 standard simulation and the RAQMS 

23 estimates of ozone export during INTEX-A, some discussion is warranted. 
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1 The main reason for the large differences between the current estimate of US 

2 photochemical ozone export during INTEX-A and the L98 standard simulations results is 

3 the anomalously cold surface temperatures during August 2004, which actually resulted 

4 in net photochemical ozone loss within the Continental US domain during the first 2 

5 weeks of August. If we restrict our budget calculations to July 01-15, we obtain a net 

6 export of photochemically produced ozone of 1.4 Tg. This export is in good agreement 

7 with the ozone export that would be obtained by applying the L98 JJA seasonal rate of 

8 1.8 Gmol/day from the standard simulation over this same period ( 1 .3Tg). However, the 

9 RAQMS ozone budget includes ozone production above the continental boundary layer 

10 while the L98 does not. The accumulated ozone P-L within the budget domain from July 

11 1-15 is 3.84 Tg. This is 85% of the accumulated continental boundary layer P-L that 

12 would be obtained for the same 15 day period using seasonally averaged P-L rates from 

13 the L98 standard simulation. The RAQMS simulation removes 2.23 Tg, or 58% of the 

14 ozone produced over the Continental US due to dry deposition during the period from 

15 July 1-15. In contrast, dry deposition removes only 33% of the ozone produced over the 

16 Continental US in the L98 standard simulation. 

17 As shown in section 3.3, the RAQMS estimates in P-L are in good agreement 

18 with observationally constrained photochemical box model estimates during INTEX-A, 

19 indicating that the current estimates of P-L are reasonable. Talbot et al., [2005] provide 

20 estimates of nocturnal ozone dry deposition during the summer based on 3 years (2001- 

2 1 2003) of ozone measurements at the Harvard Lorest site. They find median nocturnal 

22 deposition rates of 1 1 ppbv/night, which are considered representative of heavily forested 

23 regions in New England. This estimate compares very well with median RAQMS 
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1 nighttime (00Z-12Z) averaged ozone deposition velocities over New England (1 1.29 

2 ppbv/night) during INTEX-A. 

3 Taken as a whole, these comparisons indicate that the photochemical ozone 

4 export from the Continental US budget domain during July 1-15, 2004 was consistent 

5 with the L98 standard simulation results assuming that the near field ozone production 

6 due to NOy export through the continental US boundary layer was small. However, due 

7 to anomalously cold surface temperatures and resulting net ozone destruction during the 

8 first two weeks of August, the export of photochemically produced ozone was 

9 insignificant compared to the export of stratospheric ozone in the upper 

10 troposphere/lower stratosphere over the Continental US during the overall INTEX-A time 

1 1 frame (July 0 1 -August 15, 2004). 

12 The agreement between the RAQMS INTEX-A and J98 based estimates of 

13 Continental US NOy export over the INTEX-A time frame is quite good. Recall that the 

14 INTEX-A Continental US NOy budget shows accumulated NOy emissions of 0.94 Tg 

15 nitrogen and accumulated depositional loss of 0.69 Tg nitrogen (.47 Tg wet, .22 dry), 

16 resulting in a net export of 0.23 Tg of nitrogen and an export efficiency of 24%. Applying 

17 the rates from L1998 NOy budget estimates to the INTEX-A time period results in 0.86 

18 Tg of nitrogen emissions, accumulated depositional loss of 0.64 Tg nitrogen (.24 Tg wet, 

19 0.4 dry), net export of 0.23 Tg of nitrogen and an export efficiency of 27%, all of which 

20 are within 10% or less of the INTEX-A estimates. Wet deposition accounts for the 

2 1 majority of the NOy depositional loss based on the RAQMS budget calculations where as 

22 dry deposition accounts for the majority of the L1998 depositional loss. This is to be 

23 expected since the RAQMS budget domain includes the entire troposphere and 
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1 consequently the full vertical extent of wet deposition within convective cells is included 

2 in the NOy budget. The RAQMS budget analysis indicates that NOx+PAN accounts for 

3 54% of the NOy exported out of the Continential US during INTEX-A, which is 15% 

4 lower than the LI 998 estimate of 63%. This difference is due to the additional 

5 contributions from HN03 export in the middle world which is included in the RAQMS 

6 NOy budget calculations. Both INTEX-A and J1998 export efficiencies are slightly 

7 higher than the L2004 Eulerian estimates of 20% during September 1997. 

8 

9 7. Summary and Conclusions 

10 We have used aircraft, satellite, surface, and ozonesonde measurements to verify a 

11 6 week RAQMS simulation of the unified troposphere-stratosphere chemistry during the 

12 INTEX-A time period. These verification studies show that RAQMS captures the main 

13 features of the global tropospheric distribution of ozone, carbon monoxide, and NOy with 

14 reasonable fidelity, although RAQMS underestimates the median tropospheric N02 

15 distribution relative to SCIAMACHY measurements and overestimates the impact of the 

16 Alaskan wild fires on column CO relative to MOPITT. Comparisons with insitu airborne 

17 measurements shows that RAQMS reproduces the statistical characteristics of the insitu 

18 observations (median and variances) with reasonable accuracy (generally within 20%) 

19 although RAQMS tends to overestimate stratospheric influences and underestimate 

20 convective influences in the upper troposphere over the continental US (high biases in 

21 ozone and HN03 and low biases in ozone P-L and N02 above 400mb). Based on 

22 comparisons with ozonesondes from the IONS network, the assimilation of satellite based 

23 profile and column ozone measurements has been shown to have a positive impact on the 
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1 RAQMS upper tropospheric/lower stratosphere ozone analyses (mean biases of 20%), 

2 particularly during the period when higher density SAGE III limb scattering 

3 measurements were available over the Continental US. Comparisons with surface ozone 

4 measurements from the US EPA AIRNOW network show that the RAQMS surface 

5 ozone analysis captures the daily variability in surface ozone over most of the eastern US 

6 very well, with correlations between 24hr averaged measurements and the RAQMS 

7 analysis generally near 0.8. However, due to local variations in topography and 

8 emissions, the daily correlations over the Central Appalachians are considerably lower 

9 (0. 2-0.4). The RAQMS surface ozone analysis shows a systematic high bias (18ppbv at 

10 night, 15ppbv during the day) relative to AIRNow surface measurements, which is 

1 1 attributed to underestimates in nocturnal titration due to underestimates of surface NOx in 

12 urban environments. 

13 Eulerian ozone and NOy budgets during INTEX-A show that the majority of the 

14 Continental US export occurs in the upper troposphere/lower stratosphere poleward of 

15 the tropopause break. The localized ozone and NOy export was shown to occur due to 

16 convergence of tropospheric and stratospheric air in this region. These results suggest 

17 that providing a robust assessment of the influence of the Continental US on the global 

1 8 environment requires accurate representation of the long-range transport and mixing 

19 processes within this region. Continental US photochemically produced ozone was found 

20 to be a minor component of the total ozone export, which was dominated by stratospheric 

2 1 ozone that was diabatically transported into the middle world during INTEX-A. The 

22 unusually low photochemical ozone export is attributed to anomalously cold surface 

23 temperatures during the latter half of the INTEX-A mission. Efficient boundary layer 
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venting associated with cold air outbreaks during late July and mid August tended to 
reduce accumulation of ozone precusors resulting in net ozone loss during the first 2 
weeks of August. Eulerian NOy budgets during INTEX-A where shown to be very 
consistent with previously published estimates. The NOy export efficiency was estimated 
to be 24%, with NOx+PAN accounting for 54% of the total NOy export during INTEX- 
A. 
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Figure 1 : Latitude-time series of the frequency of observations used in the RAQMS 
assimilation. Symbols indicate the location of solar occultation and limb scattering 
observations. Contours indicate the density (zonal mean % of total) of cloud-cleared total 
column measurements. 
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Figure 2: Comparison between RAQMS and cloud cleared MOPITT CO column (10 18 
mol/cm 2 ) for the period from July 01-August 15, 2004. Continental US budget domain is 
indicated in white. 
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Figure 3: Comparison between RAQMS and cloud-cleared SCIAMACHY tropospheric 
N02 column (10 15 mol/cm 2 ) for the period from July 01-August 15, 2004. Continental 
US budget domain is indicated in white. 
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Figure 4: Comparison between RAQMS and climatological (1979-2000) tropospheric 
ozone column (DU) during the period from July 01- August 15, 2004. Continental US 
budget domain is indicated in white. The location of the WMO ozonesonde station at 
Santa Cruz, Tenerife is indicated by a diamond. 
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Figure 7: Comparison of time averaged IONS ozonesonde and coincident RAQMS ozone 
analyses during July 01 -August 15, 2004. 
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Figure 8: Comparisons between RAQMS and EPA AIRNOW surface ozone 
5 measurements during July 01 -August 15, 2004. 
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Figure 9: Comparisons between RAQMS and INTEX-A DC8 insitu observations of CO, 
03 (upper panels), N02, Total PANs (middle panels), HN03, and observationally 
constrained boxmodel P-L calculations during July 01 -August 15, 2004. CO and 03 are 
in ppbv, N02, Total PANs, and HN03 are in pptv, P-L is in ppbv/day. 


56 


1 

2 

3 



6 

7 

8 Figure 10: Zonally averaged global 380K (upper panels) and cross tropopause (lower 

9 panels) fluxes of 03 (10 11 mol/cm 2 /sec, left panels) and NOy (10 8 mol/cm 2 /sec, right 

10 panels) during July 01 -August 15, 2004. The thin solid line is the net flux, which is the 

1 1 sum of isentropic (dash-dot), diabatic (dotted), and movement of the surface (dashed). 

12 The bold line denotes the mean altitude (km) of the 380K and tropopause surface. 
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Figure 1 1 : Zonally averaged Continental US 380K (upper panels) and cross tropopause 
(lower panels) fluxes of 03 (10 11 mol/cm 2 /sec, left panels) and NOy (10 8 mol/cm 2 /sec, 
right panels) during July 01 -August 15, 2004. The thin solid line is the net flux, which is 
the sum of isentropic (dash-dot), diabatic (dotted), and movement of the surface (dashed) 
The bold line denotes the mean altitude (km) of the 380K and tropopause surface. 
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Figure 12: Time averaged zonal mean distribution of Continental US ozone (ppbv), Net 
P-L (ppbv/day), convective mixing tendencies (ppbv/day) and absolute assimilation 
increment (%) as a function of altitude (km) and latitude within the Continental US 
budget domain during July 01 -August 15, 1004. 
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Figure 13: Time averaged ozone number densities (mol/cm3) for each of the lateral 
boundaries of the Continental US budget domain for July 01-August 15, 2004. 
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Figure 14: Time averaged ozone fluxes (mol/cm2/sec) for each of the lateral boundaries 
of the Continental US budget domain for July 01 - August 15, 2004. Negative fluxes are 
into the domain. 


61 



1 

2 

3 

4 

5 


Accumuloted 03 July 01 -Aug 15, 2004 (CONUS SFC-380K) 



Julian Doy 

6 

7 

8 

9 Figure 15: Time series of accumulated changes in Continental US ozone (Tg) for July 01 

10 - August 15, 2004 due to P-L, 380K fluxes, lateral fluxes, ozone assimilation, and dry 

1 1 deposition. The actual and computed accumulation is also shown. 
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Figure 16: Time averaged zonal mean distribution of Continental US NOy (ppbv), Net 
production due to lightning NOx (ppbv/day), convective mixing tendencies (ppbv/day) 
and wet deposition (ppbv/day) as a function of altitude (km) and latitude within the 
Continental US budget domain during July 01 -August 15, 1004. 
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Figure 17: Time averaged NOy number densities (mol/cm3) for each of the lateral 
boundaries of the Continental US budget domain for July 01-August 15, 2004. 
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Figure 18: Time averaged ozone fluxes (mol/cm2/sec) for each of the lateral boundaries 
of the Continental US budget domain for July 01 - August 15, 2004. Negative fluxes are 
into the domain. 
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7 Figure 19: Time series of accumulated changes in Continental US NOy (Tg Nitrogen) for 

8 July 01 - August 15, 2004 due to emissions, 380K fluxes, lateral fluxes, lightning NOx, 

9 wet and dry deposition. The actual and computed accumulation is also shown. 
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